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O '■ Abstract 

a: 

O 1 We present a general method for calculating inclusive cross sections for 

i-^ ■ the production of heavy quarkonium states with definite polarization within 



the NRQCD factorization approach. Cross sections for polarized production 



5— i 

■ can involve additional matrix elements that do not contribute to cross sec- 

tions for unpolarized production. They can also include interference terms 
between parton processes that produce QQ pairs with different total angular 
momentum. The interference terms cancel upon summing over polarizations. 
Our method can be generalized to N dimensions and is therefore compatible 
with the use of dimensional regularization to calculate radiative corrections. 
We illustrate the method by applying it to the production of J/tp via the 
parton processes qq — ► cc and gg — ► cc. 
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I. INTRODUCTION 



Calculations of inclusive production rates of heavy quarkonium have recently been placed 
on a firm theoretical foundation by the development of a factorization approach based on 
nonrelativistic quantum chromodynamics (NRQCD) In this formalism, the cross sec- 
tion is expressed as a sum of products of short-distance coefficients and NRQCD matrix 
elements. The short-distance coefficients can be calculated as perturbation series in the 
coupling constant a s at the scale of the heavy quark mass. The matrix elements scale in a 
definite way with the typical relative velocity v of the heavy quark in the quarkonium state. 
Thus the production cross sections can be calculated systematically to any desired order in 
a s and v 2 in terms of well-defined NRQCD matrix elements. There have been many recent 
applications of this formalism to quarkonium production in various high energy processes 

Most calculations of quarkonium production have been carried out using a covariant pro- 
jection formalism developed for calculations in the color-singlet model |3[. In the amplitude 
for producing a cc pair with total momentum P and relative momentum q, the spinor factor 
uv is replaced by an appropriate Dirac matrix that projects out a state with the desired 
angular-momentum quantum numbers. The resulting prescription is relatively simple for S- 
wave states. The Dirac matrix is proportional to 75 — 2m c ) for a 1 So state and tfffi — 2m c ) 
for a 3 S\ state, and the relative momentum q is set to zero. For P-waves, the Dirac matrix 
is more complicated, and the amplitude must be differentiated with respect to the relative 
momentum. The resulting expression for the amplitude can be very complicated. 

The covariant projection method also has a number of other drawbacks. For each angular 
momentum state, a separate calculation of the cross section is required, beginning at the 
level of the amplitude for cc production. Another drawback is that it is difficult to generalize 
the projection method to allow the calculation of relativistic corrections to the cross section. 
There is also a potential difficulty in calculating radiative corrections using this method. 
The Dirac matrices that are used to project onto the appropriate angular momentum states 



2 



are specific to 3 space dimensions. However the most convenient method for regularizing 
the ultraviolet and infrared divergences that arise in higher order calculations is dimensional 
regularization. It is not easy to generalize the projection matrices to N dimensions, since 
the representations of the rotational symmetry group are different for each integer value of 
N. 

In this paper, we develop a new method for calculating quarkonium production rates 
that fully exploits the NRQCD factorization framework. The short-distance coefficients are 
calculated in a form that holds for every quarkonium state. The corresponding NRQCD 
matrix elements are then simplified using rotational symmetry and the approximate heavy- 
quark spin symmetry of NRQCD. It is only at this stage that the angular momentum 
quantum numbers of the quarkonium state come into play. Relativistic corrections are easily 
calculated in this framework. Our method for calculating the short-distance coefficients is 
readily generalized to N spacial dimensions, so that dimensional regularization can be used 
to cut off infrared and ultraviolet divergences. After removing the divergences, one can 
specialize to N = 3 and use rotational symmetry and spin symmetry to simplify the NRQCD 
matrix elements. This approach allows the consistent use of dimensional regularization to 
calculate inclusive heavy quarkonium production and decays. 

Our method has important implications for the production of quarkonium states with 
definite polarization. We demonstrate that cross sections for polarized production involve 
new matrix elements that do not appear in cross sections for unpolarized production. Thus, 
measurements of cross sections for unpolarized production are not necessarily sufficient to 
predict production rates for polarized quarkonium states. We also show that cross sections 
for polarized production can involve interference between parton processes that produce 
cc pairs with different total angular momentum J. The interference terms vanish upon 
summing over polarizations. Thus the cross sections for producing cc states with definite 
total angular momentum are not sufficient to determine all the short-distance coefficients in 
cross sections for polarized production. 

The NRQCD factorization approach is summarized in section II. In section III, we present 
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a general matching procedure for calculating the short- distance coefficients in the factoriza- 
tion formula. We illustrate the method by applying it to the parton processes qq — > cc and 
gg — > cc. The resulting expressions for the cross sections hold for every quarkonium state. In 
section IV, we show how the NRQCD matrix elements can be simplified by using rotational 
symmetry and the approximate heavy-quark spin symmetry of NRQCD. For simplicity, we 
focus on the matrix elements for J /if) production. In section V, we calculate the cross section 
and the spin alignment of the ip from the parton processes qq — > cc and gg — > cc. We point 
out that several previous calculations of quarkonium production need to be reconsidered in 
the light of our results. 



II. NRQCD FACTORIZATION FORMALISM 

We consider the inclusive production of a quarkonium state H with momentum P and 
helicity A via a parton process of the form 12 — > H(P, A) + X. The differential cross section, 
summed over additional final states X, can be written 



]TAx(12 -> H(P,X)+X) 
x 

1 d 3 P 
~ AE x E 2 v 12 (2jc) 3 2E p v 



J2(^) 4 6 4 (h + k 2 -P-k x ) |T 1M(PiA)+x | 2 , (2.1) 



where Ep = J M\ + P 2 and the sum over X includes integration over the Lorentz-invariant 
phase space for the additional particles. 

The NRQCD factorization formalism can be used to factor the cross section ( p.l|) into 
short-distance coefficients and long-distance matrix elements [0: 

1 ii 3 P 

$>a(12 - A) + X) = -— — Y.C mn {PMM){0^). (2.2) 

x ^E x E 2 v X 2 {2-k) 6 2E p mn 

The coefficients C mn are functions of the kinematic variables P, ki, and k 2 . They take 
into account the effects of short distances of order l/m c or smaller, and therefore can be 
calculated as perturbation series in the QCD coupling constant a s (m c ). The matrix elements 
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(Omn^) are expectation values in the NRQCD vacuum of local 4-fermion operators that have 
the structure 

= x^'liJ V H{ x) tfKnX, (2-3) 

where tp and \ are the field operators for the heavy quark and antiquark in NRQCD, and K, n 
and K. m are products of a color matrix (1 or T a ), a spin matrix (1 or a 1 ), and a polynomial 
in the gauge covariant derivative D. The projection operator Vm\) can be written 

V HW = ^2 \H(P = 0, A) + S){H(P = 0,X) + S\. (2.4) 

s 

The sum is over soft hadron states S whose total energy is less than the ultraviolet cutoff A 
of NRQCD. Thus this operator projects onto the subspace of states which in the asymptotic 
future include the quarkonium state H(X) at rest plus soft hadrons. The normalization of 
the meson states in (|2.4j) must coincide with those in the T-matrix element in (2.1). The 
standard relativistic normalization is 

(H(P',X')\H(P,X)) = 2E P (2nf6 3 (P-P')6 xy . (2.5) 

With this normalization of states, the projection operator Vh(x) has energy dimension -2. 
If the colliding particles are leptons, the cross section for quarkonium production is given 



directly by the factorization formula (|2.2j ). If the colliding particles are hadrons, the cross 
section ( |2.2|) must be folded with parton distributions for the partons 1 and 2 in the colliding 
hadrons. In this case, the derivation of the factorization formula requires that the transverse 
momentum p? of the quarkonium be large compared to Aqcd, the scale of nonperturbative 
effects in QCD. This restriction to large pt follows from the diagrammatic analysis that 
underlies the factorization formula . This analysis shows that the dominant contributions 
to the cross section can be factored into (a) hard-scattering amplitudes for parton processes 
of the form 12 — > cc + 34 . . ., (b) jet-like subdiagrams for the incoming partons 1 and 2 and 
the outgoing partons 3, 4, etc., (c) a subdiagram involving a cc pair with relative momentum 
that is small compared to the quark mass m c , and (d) a soft part. The soft part involves soft 



gluons that couple to the jet-like subdiagrams and to the cc subdiagram, but not to the hard- 
scattering subdiagram. The effects of soft partons that are exchanged between the various 
subdiagrams cancel upon summing over all possible connections of the soft partons. This 
cancellation is effective provided that the cc pair has large transverse momentum relative 
to the incoming hadrons. All the effects of soft partons can then be factored into parton 
distributions associated with the incoming hadrons, fragmentation functions associated with 
the outgoing partons produced by the hard scattering, and a factor associated with the cc 
pair that depends on their relative momentum q. In this step in the derivation of the 
factorization formula, the short-distance scales pt and m c are separated from the long- 
distance scale Aqcd- The effects of the short-distance scales appear only in the hard- 
scattering amplitudes, and all effects of the scale Aqcd are factored into parton distributions, 
fragmentation functions, and cc factors. 

The remaining step in the derivation of the factorization formula involves separating the 
scale m c in the hard-scattering amplitude from the scale m c v of the relative momentum 
in the charmonium state. This separation can be accomplished by Taylor-expanding the 
hard-scattering amplitude in powers of q. A naive Taylor expansion generates ultraviolet 
divergences in the cc factor. The effective field theory framework provided by NRQCD 
makes this Taylor expansion meaningful. The long-distance factors generated by the Tay- 
lor expansion can be identified with matrix elements of local operators in NRQCD. The 
renormalization framework of NRQCD allows the ultraviolet divergences in these matrix 
elements to be systematically removed. The resulting expression for the cross section has 



the form of the NRQCD factorization formula ( |2.2| ). The effects of the scale m c appear only 
in the short-distance coefficients C nm . All effects involving momentum scales of order m c v 
or smaller have been factored into the NRQCD matrix elements {O^j^). 

The breakdown of the factorization formula at small when the colliding particles are 
hadrons can be attributed to "higher- twist" processes. In the factorization formula ( |2.2| ), the 
c and c that form charmonium are assumed to be produced by a hard-scattering collision 
involving a single parton from each of the colliding hadrons. Higher-twist processes can 
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involve more than one parton from a single hadron. An example of a higher-twist process 
for charmonium production is one in which the c and c that bind to form charmonium are 
produced by a quantum fluctuation of the hadron into a state containing a cc pair [|]]. The 
transition of the cc to a charmonium state can then be induced by the collision with the 
other hadron. The short- distance part of this process is the fluctuation of the hadron into a 
state containing the cc pair. The cross section for the scattering of the cc pair from the other 
hadron involves long-distance effects that cannot be expressed in terms of matrix elements 
of local NRQCD operators. 

III. SHORT-DISTANCE COEFFICIENTS 

The short-distance coefficients C mn in (|2.2|) can be determined by matching perturbative 
calculations in full QCD with the corresponding perturbative calculations in NRQCD. In 
most previous calculations, the coefficients were determined by matching cross sections for 
producing cc pairs with definite total angular momentum. However, cc states with different 
angular-momentum quantum numbers can have a nonzero overlap with the same final state 
\H + S). The resulting interference terms cannot always be obtained by matching cross 
sections. Below we present a matching prescription which is sufficiently general to provide 
these interference terms. 

A. General matching prescription 

Let cc(q, £, rj) represent a state that consists of a c and a c that have total momentum P, 
3-momenta ±q in the cc rest frame, and spin and color states specified by the spinors £ and 
rj. These spinors are 2-component Pauli spinors with a color index. Using the abbreviated 
notation cc = cc(q, £, rf) and cc' = cc(q', 77'), the matching condition is 

]T(27r) + k 2 - p - k x ) V 2 _ c t +X r 12 _ cc - +x 

X 

= E C mn(P, h, k 2 ) (X^'LV ^,c, ^IC nX )- (3.1) 

mn 
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The operator V C c',cc in the matrix element in (|3.1|) is defined by 

Vcc',cc = Y,\ c (4,£)c(-<l!,v') + S)(c(q,0c(-<i,v)+S\. (3.2) 

The sum is over soft parton states whose total energy is less than the ultraviolet cutoff A 
of NRQCD. The normalization of the cc states in Q3.2| ) must coincide with those in the 
T- matrix elements in fl2.1|) . The standard relativistic normalization is 

(c(qU0%WMqi,£)c(q 2 ,77)> = 4E qi E q2 (27r) 6 5 3 ( qi - q' 1 )5 3 (q 2 - q' 2 ) ^ (3.3) 

where the spinors are normalized so that = rfrj = 1, and similarly for £' and rf . In 
expressions like both the spin and color indices are contracted. Note that with the 
normalization ( |3 . 3| ) , the operator P cc / iCC in ( |3.2| ) has dimension -4. The difference in the 
dimensions of the operators P# and P cc ' jCC matches the difference in the dimensions of the 
T- matrix elements in Q2.ip and ( |3.1| ). 

To carry out the matching procedure, the left side of (|3.1[ ) is calculated using perturbation 
theory in full QCD, and then expanded as a Taylor series in q and q'. The matrix elements 
on the right side of ( |3.1|) are calculated using perturbation theory in NRQCD, and then 
expanded as Taylor series in q and q'. The short- distance coefficients C mn are obtained by 
matching the terms in these Taylor expansions order by order in a s . 

B. Matching for qq — ► cc 

We first illustrate the matching procedure by applying it to the process qq — > cc. The 
T-matrix element in Feynman gauge is 

T 1M = g 2 ^ v(k 2 ) llM T a u(h) u(p)^T a v(p). (3.4) 

Dirac and color indices are implicit in the spinors. The 4-momenta of the outgoing c and c 
can be expressed as p = ~P + Lq and p = ~P — Lq, where P is the total 4-momentum of the 
cc pair, q is the relative 3-momentum of the c in the center-of-momentum (CM) frame of the 
pair, and V Jj i is a Lorentz boost matrix. The boost matrix V Jj i transforms a purely spacelike 
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4- vector, such as (0, q), from the CM frame where the components of P are (2E q , 0) to the 
frame in which its components are P M . The matrix is given explicitly in Appendix A. The 
expressions for the Dirac spinors u(p) and v (p) in terms of Pauli spinors £ and rj are also 
given in Appendix A. Identities in Appendix A can be used to expand u(p) , y fl T a v(p) as a 
Taylor series in q. To linear order in q, the T-matrix element ( J3.4J) is 

T 1M = -jf- vikzh^uih) L\ £VT> (3.5) 

The components of the boost matrix to linear order in q are 

L\ = -L/-, (3.6a) 

L\ = S ji -P j P i + — PP*, (3.6b) 

2m r 



where Pp = y4 m c + P 2 . 

In order to carry out the matching procedure, we first calculate the left side of ( |3.1|) . In 
this case, there is no sum over X. Multiplying ( |3.5| ) by the complex conjugate of T^^cc 1 
and averaging over initial spins and colors, we obtain 



(2vr) 4 5 4 (A; 1 + k 2 - p) £ r x 



12-+C5 



A 2 2 

= (2tt) 4 5 4 (A; 1 + A; 2 - P) [<P - n^'n*] ^^T^^a^r). (3.7) 

Q L J 

We have simplified the expression by using the identities (|A5|) and (|A6aj) in Appendix A. 



We have also used the fact that (k\ — k'i) v JJ l i = —VP^n 1 , where n is a unit vector, which 
follows from ( |A6b ). In a frame in which P, k 1; and k 2 are collinear, they are also collinear 



with n. 

We next consider the right side of the matching equation ( |3.1|) . We must construct an 



operator of the form (|2.3|) whose vacuum matrix element, when calculated to leading order 
in a s , reproduces the spinor factor in ( |3.7| ) upon expanding to linear order in q and q'. One 
can see by inspection that the appropriate operator is x^ j T a ^V C c',cc^^T a x- At leading 
order in a s , only the cc states in the projection operator P cc / jCC contribute: 

(xVT>P cc -, cc -V>VT a X > = (OlxV'T^lc^e'^')) (cc(q,tv)\^^T a x\0)- (3.8) 



The vacuum-to-cc matrix element reduces at leading order in a s to 



(cc(q,e,r/)|^VT a X |0) = 2E q £a i T% 



(3.9) 



where the factor of 2E q arises from the relativistic normalization of states. Expanding to 
linear order in q and q', the matrix element ( p.8| ) reduces to 



(xVT> P cc - Vc - ^a'T a X ) = 4m c 2 t/Vt^'CVT^ 



(3.10) 



Comparing (|3.7|) and ( |3.10|) , we can read off the short-distance coefficient for the matrix 
element: 



C Jl = (2vr)M 4 (fc 1 + k 2 - P) 



9m 2 . 



5 ji - h j ti 



(3.11) 



Having determined the short- distance coefficient, we can insert it into ( |2.2|) to obtain an 
expression for the differential cross section for producing a meson H with helicity A via the 
parton process qq — > cc: 

1 d 3 P 



da(qq^H(P,X)) 



AE 1 E 2 v 12 (2nf2Ep 



(2vr) 4 5 4 (A; 1 + k 2 - P) 



x 



9m 2 



5 ji - n j ri] (x^ j T a 4> V H (X) ^Vr°x). (3.12) 



Using the delta function to integrate over the phase space of P, we find that the cross section 



is 



a(qq ^ H(\)) = S(s - 4m; 



2^ ^ol 2 s 



36m 4 



( X *a j T a iP V H (X) VVr°x>, (3.13) 



where s = (&i + /C2) 2 and we have assumed that the 3-momentum of the colliding partons 
is along the z-axis. Note that the expression ( |3.13| ) is dimensionally correct. Each of the 
quark fields has dimension | and the projection operator has dimension -2, so the dimension 



of the matrix element is 4. The expression ( |3.13| ) gives a contribution to the cross section 
for every quarkonium state H. 
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C. Matching for gg — ► cc 



We next illustrate the matching procedure by applying it to the process gg 
T-matrix element can be decomposed into 3 independent color structures: 



T 



V 6 2 2 



cc. The 



(3.14) 



The Dirac factor for the 5 term is 



= u{p) 



v(p). 



(3.15) 



The term D^ vc in (|3.14j) differs from (|3.15|) only by inserting the color matrix T c between 
the spinors. The Dirac factor for the f abc term in Feynman gauge is 



F^ c = u(p) 



2p ■ k\ 2p ■ k 2 

o (P + k 2 YY ~(P + hTr + g» v (fa- fa) 



p2 



T c v{p). (3.16) 



Using the identities in Appendix A, we expand ( |3.15|) and ( |3.16| ) to linear order in q: 



2mi 



e^ x »(h - k 2 ) x P p g v 



+ 



-^(h ■ L) n ((h - k 2 yU 3 + (h - k 2 YL^) - \{k, ■ L),{P^L\ - P»V n ) 



mi m r K 



q n ^ j V, (3.17) 



F^ c = ~ (k^L v j - k»L^) ^a j T a V - ^ e^ xp (h - k 2 ) x P p (h ■ L) n q n (3.18) 

The first term in F^ vanishes when contracted with the polarization vectors e p (ki) and 
e b u (k 2 ) in fl3.14| ). Multiplying (|3.14j) by the complex conjugate of T gg ^ cE ' and averaging over 
the spins and colors of the initial gluons, we obtain 



ET* T 
1 gg^cc l £ g 



+ 



m 



(S mn - fTh n )(6 ij - hW) + (S mj - h m n j )(S ni - fFh 1 ) 
+ S mi S nj _ (S mi - n m h%8 nj - n n h j ) 
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27 1 

+ n m n n g"V ri^T a £^T a ri }. (3.19) 

8m^ J 

At leading order in a s and to linear order in q and q', the spinor factors on the right side 
of ( |3.19|) can be identified with the following matrix elements: 



(xVTW-^x) = Aml^^V, (3.20a) 

( X j T a ^V c5 ,, c - c ^T a x) = Am 2 c ^Tr^T% (3.20b) 

( x \-0™)T a ^V^, cB ^{-0 n )T a X ) = Am 2 c q' m q n ^T a C^T a V , (3.20c) 

(x\-0 m )a i T a i)V C c>, C c^\-iD n )a j T a x) = Am 2 c q' m q n ^o l T a £,'£}o j T a ri, (3.20d) 

where ip^~Dx = i'^X ~ (D^/O^X- 

Matching terms on the left side and right side of ( p.l| ), we determine the short- distance 
coefficient C mn for each of the matrix elements. Inserting these coefficients into ( |2.2| ), we 
obtain the differential cross section for production of H(X) from the parton process gg — > cc. 
The final result, after integrating over phase space, is 

7r 3 (T 2 f 1 ^ 

a{gg^H{\)) = 8(s - Am 2 c )-^UxHV H(x) tf X ) + - (x^V H(x) ^T a x ) 



+ — 



(S mn - z m z n )(S ij - z l z j ) + (S mj - z m z J )(5 ni - z n z l ) + 5 mi S nj 

_ ( 6 mi _ $m$i)(Sn3 _ ^i)] (( % t (_ m )a^ V H (X) ^{-0 n )o j x) 

1 ^ <-» <-» 

o 

27 <-+ +-» 1 

+ — ( x t ( _| Dm)T a^p H(A) V;t( _| jD n )T a x) (g^) 



The expression ( |3.21| ) for the cross section applies equally well to any quarkonium state H. 
The relative importance of the various terms for a given state H depends on the magnitude 
of the matrix elements. 
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IV. REDUCING THE MATRIX ELEMENTS 



The cross sections (|3.13| ) and ( p.21|) hold for any quarkonium state H, with the depen- 
dence on H appearing only in the matrix elements. The matrix elements can be simplified 
by using rotational symmetry and the approximate heavy-quark spin symmetry of NRQCD. 
Their magnitudes can be estimated using velocity-scaling rules for NRQCD matrix elements. 
In this section, we apply these methods to simplify the matrix elements that appear in in- 
clusive cross sections for the production of J/ip. The extension to other quarkonium states 
will be presented elsewhere ||. 



A. Rotational symmetry 



Under rotations, each of the matrix elements in ( [3.13 ) and ( |3.21 ) transforms as a com 



ponent of a tensor. The indices of the tensor include vector indices, such as i, j, m, and n, 
and the helicity A in the projection operator Vh(x)- The helicity A appears both in the ket 



and in the bra in the expression for the projection operator given in (2^). If the meson has 
total angular momentum J, then Vh(\) transforms like the (A, —A) component of the direct 
product representation J_® J_. Thus the entire operator transforms as an element of the 
representation J_® J_® 1® 1® where there is a 1 for each vector index of the operator. 
The vacuum matrix element of the operator can be expressed in terms of iV independent 
matrix elements, where N is the number of times the trivial representation appears in the 
decomposition of the direct product J® J_® 1 <S> 1® ... into irreducible representations. 

In the case of the J/if>, the independent matrix elements can be determined by elementary 
tensor analysis. The ijj has total angular momentum J — 1, and the helicity label A is just 
a vector index in a spherical basis. The spherical basis and the Cartesian basis are related 
by a unitary transformation with matrix 
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u 



Xi 



^ -1/V2 -i/y/2 ^ 
1 

1/V2 -i/y/2 



(4.1) 



Xi 



The helicity labels A in the ket and in the bra in (|2.4j ) are related to the corresponding Carte- 
sian indices % and j by the transformation matrices Ux% and Uj x , respectively. Since there 
are no momentum vectors that the matrix element can depend on, it must be expressible in 
terms of the invariant tensors 5 mn and e lmn . The number of independent matrix elements is 
the number of invariant tensors that can be formed out of the available indices. 

We first consider operators with no vector indices, such as (x ^"Pfix)^ x) ■ The only 
indices that are available are A and A from the projection operator Vmx)- The only invariant 
tensor is 5\\ = 1. Therefore, we have 

{xHV m ^x) = 3<xtyP^ + X>, (4.2) 

where V$ = Sa^Vo)- Standard production matrix elements Of( 2S+1 Lj) and 0$ ( 2S+1 Lj) 
were defined in Ref. [|l] using a projection operator analogous to Q2.4p , but constructed out 
of states with the standard nonrelativistic normalization. As discussed in Appendix B, the 
two projection operators differ at leading order simply by a normalization factor of 4m c . 
Thus, we can write 

(xV^(A) ^X) = \rn c (OfCSo)), (4.3a) 
&T^V^x)^T a x) = ^m c (OtCS )), (4.3b) 

where the matrix elements on the right sides are defined in Appendix B. 

We next consider matrix elements with two vector indices, such as (x ^V^Vca)^ a " 'x) ■ 
Rotational symmetry implies that such a matrix element can be expressed in terms of 3 
independent tensors: b l \ U\iUj x , and U\jU\ x . Thus there are three independent matrix 
elements. If we sum over the helicities A, the only possible tensor is , so there is only one 
possible matrix element. Thus, up to corrections of relative order v 2 , we have 
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£<XVV7V (A) V-Vx> = I ^ m c (Of C 3 ^)), (4.4a) 
^(xV^^A) ^V'T-x) = ^ 5ij m < (°tCSi)), (4.4b) 

A 6 

^(-0^V m ^(-^) X ) = U^m c (Of( 1 P 1 )), (4.4c) 

A 6 

^{X ] {-{D i )T a ^V^ {X )^ {-{D j )T a x) = U^m^OfCPx)), (4.4d) 

A 6 

where the matrix elements on the right sides are defined in Appendix B. In the cross section 
( p.21j ), we also require the matrix element {x\—\D^)T a ipVf{\)ip\—\D^) T a x) contracted 



with z x z? . This can be expressed in terms of (Of { l Pi}) and one additional matrix element: 

35ao - 1 



+ <xK-{&)T a $V m il>\-±I?)T*x). (4-5) 



Upon summing over helicities, we recover ( |4.4d| ) 



Matrix elements with 4 vector indices, such as (x^cr t (—^D m )ipT > H(x)^o' : '(—^D n )x), can 
be reduced by rotational symmetry to 15 independent matrix elements. Upon summing over 
helicities, the number of independent matrix elements is reduced to 3. 



B. Vacuum-saturation approximation 

The vacuum matrix element (x o'V'^VmV' ^x) can De simplified by using the vacuum- 
saturation approximation. In this approximation, only the NRQCD vacuum is retained in 
the sum over soft states in the projection operator ( |2.4j) . The matrix element then reduces 
to 

(xVvp^Vx) « (o|xVv|V(p = o,a))(v(p = o,a)|vVx|o). (4.6) 

As pointed out in ref . |l| , the vacuum-saturation approximation is a controlled approximation 
in the case of heavy quarkonium, with a relative error of order v 4 . The reason for this is that 
the matrix element involves a transition from a color-singlet cc state to states of the form 
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\ip{\) + S), where S is a soft state. Since the if) is a color singlet, the state S must also be 
a color singlet. Suppose it is not the NRQCD vacuum. Since a single chromoelectric dipole 
transition will produce a color-octet state, at least two chromoelectric dipole transitions are 
required to produce the state S. The amplitude is suppressed by a power of v for each such 
transition. There is another power of v 2 from the complex conjugate amplitude, and this 
gives an overall suppression factor of v 4 . 

After using the vacuum-saturation approximation, the matrix element ( |4.6D can be fur- 
ther simplified by using rotational symmetry. The vacuum-to-^ matrix element must have 
the form 

(V>(P = 0,A)|^Vx|0) = ul^m^^R^. (4.7) 

The factor of ^2M^ comes from the relativistic normalization of the charmonium state. The 
factor ^j3/2n has been inserted so that can be interpreted as the nonrelativistic radial 
wavefunction evaluated at the origin. Inserting ( |4.7| ) into ( |4.6| ), it reduces to 

(xVV^(a)VVx) = U Xi Ul -M^l 2 + 0(v 4 m c \R^\ 2 ). (4.8) 



C. Heavy-quark spin symmetry 

The most powerful tool at our disposal for simplifying the production matrix elements 
of heavy quarkonium is the approximate heavy-quark spin symmetry of NRQCD. The spin 
symmetry transformations of the heavy-quark and antiquark fields are 

V»(x,t) -> V V»CM), (4.9a) 
X {x,t)^W X (x,t), (4.9b) 

where V and W are SU(2) matrices. Under a rotation, the transformation of the fields is 

ip{x,t) -> V ^(0.x,t), (4.10a) 
x(x,t)->Vx(0.x,t), (4.10b) 
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where O is the 0(3) matrix whose elements are 

O ij = -tr(yVW), (4.11) 

While rotational symmetry is an exact symmetry of NRQCD, spin symmetry is only an 
approximate symmetry. NRQCD is equivalent to full QCD in the sense that the coefficients 
of the terms in the NRQCD lagrangian can be tuned so that the effective theory reproduces 
the quarkonium spectrum and low energy quarkonium matrix elements to any desired order 
in v. Minimal NRQCD reproduces the spectrum and matrix elements to a relative accuracy 
of order v 2 . This theory has exact heavy-quark spin symmetry. With the t> 2 -improved 
NRQCD lagrangian, the accuracy is improved to v 4 . However the v 2 -improved lagrangian 
includes a term f/^B • crip — x'B • a\ that breaks the spin symmetry. Thus spin symmetry 
gives relations between matrix elements that are accurate up to corrections of relative order 
v 2 . 

The consequences of spin symmetry are particularly simple for spin-triplet matrix ele- 
ments of the ip, such as (x^ ^ipV^x)^ cr-' x) ■ Since the ip is an S-wave state, its helicity A is 
a spin index. Spin symmetry implies that the helicity label A in the ket of the projection 
operator Vmx) must be matched with the vector index i, while the A in the bra must be 
matched with the index j. Thus the matrix element must be proportional to U\iUj X . The 
proportionality constant is obtained by summing over A and contracting with 8 l K We find 
that, up to corrections of relative order v 2 , 

(xVV%) = \ U Xl U} x m c (Ot( s S l )}. (4.12a) 

( X VT>P^ (A) ^VT a X ) = ^U Xl U] x m c (Ot( 3 S 1 )}, (4.12b) 



The relative error of v 2 in ([4.12a) is larger than the relative error of v 4 obtained by using 



the vacuum saturation approximation ( |4.8| ). One can easily show that the error in ( |4.12a|) 
can be improved to v 4 by replacing 4m c by 2M^. 

The matrix element (x t (T i (-|^ ,m )^(A)^ t ^(-fD n )x) can also be reduced to a single 
matrix element by using spin symmetry. This symmetry implies that the tensor structure 
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in the indices A, i, and j must be the same as in ( f4.12b|) . Rotational symmetry then implies 
that the matrix element must also be proportional to 5 mn . Thus, up to corrections of relative 
order v 2 , we have 

l^^iD^V^ V>V(-0")x> = 4 U Xi U] x 5™ m c (Of( 3 P )), (4.13a) 
(xV J (-fB m )T>^ (A) ^V(-|£")T a X > = 4 U Xi Uj x 5 mn m c (Ot( 3 P )}. (4.13b) 

Spin symmetry can also be applied to spin-singlet matrix elements of the if>, such as 
(xK~% D % )il)V^(\)il)\— | Di)x)- A spin-singlet operator like ^(—-D^x creates a cc pair in 
a spin-singlet state. At leading order in v 2 , the transition to a final state of the form \if) + S) 
must involve the term if)<H ■ crip — x^B • ax from the v 2 -improved NRQCD lagrangian. The 
resulting constraints on the matrix element can be deduced using angular momentum theory. 
These constraints are simply those that are already provided by rotational invariance. Spin 
symmetry does give new relations between the production matrix elements of the if> and 
those of the r] c , but they are not of great practical significance. 



D. Velocity-scaling rules 

The relative importance of the NRQCD matrix elements that appear in cross sections 
for if} production are determined by how they scale with the relative velocity v of the charm 
quark in the meson |]]. When applied to the matrix elements (O n ( 2S+1 Lj)), the scaling rules 
are fairly simple. The operator O n ( +1 Lj) creates and annihilates a pointlike cc pair in the 
color state n and in the angular momentum state 2S+l Lj. The scaling of the matrix element 
is determined by the orbital-angular-momentum quantum number L and by the number of 
chromoelectric and chromomagnetic dipole transitions that are required for the cc pair to 
reach the state cc(n, 2S+l Lj) from the dominant Fock state of the meson. If the minimum 
number of these transitions is E and M, respectively, the matrix element (O n ( 2S+l Lj)) 
scales as v ^l+2E + am _ 

The dominant Fock state of the if> consists of a cc pair in a color-singlet 3 Si state. 
Thus the matrix element (Oi( 3 S\)) scales as v 3 . The cc pair can reach the color-octet 
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states cc(8, 3 Pj), cc(8, 3 Si), and cc(8, 1 5'o)) through a single chromoelectric, a double chro- 
moelectric, and a single chromomagnetic dipole transition, respectively. Thus the matrix 
elements (Os( 3 Pj)), (Os( 3 Si)), and (Og^So)) all scale as v 7 . All the other matrix elements 
(O n { 2S+1 Lj)) scale as v 11 or smaller. Those that scale like v 11 include (O^Pt)), (Og^Px)), 
(O^Sq)), and (O^Pj)). The matrix element (x j (-0 3 )T a °ipV m ^(-f D 3 )T a X ) also 
scales as w 11 like (C 8 ( 1 P 1 )). 



In this section, we combine the results of Sections III and IV to obtain the cross sections 
for production of ip with definite helicity via the order-a^ parton processes qq — > cc and 
gg — > cc. These cross sections can be folded with parton distributions to obtain cross 
sections for ip production in hadron-hadron scattering. The resulting cross sections should 
not be taken too seriously, since these parton processes produce a ip with zero transverse 
momentum. The factorization formula ( |2.2| ) is therefore not strictly applicable, since its 
derivation requires the transverse momentum of the ip to be large compared to Aqcd- We 
choose to ignore this difficulty, since our primary purpose is to illustrate the calculation of 
cross sections for polarized production. 

The cross section for xjj production via the parton process qq —>■ cc is obtained by inserting 
the matrix element ( |4.12b| ) into (|3.13|) : 



We have used the fact that the matrix U in ( [4.1|) satisfies J2i U\iU iX = 1 and J2j Uj X z^ = 5\ . 
The cross section ( |5.1|) implies that ?/>'s produced by this process are transversely polarized, 
which could be anticipated from the fact that the gluon interaction with the light quark 
conserves helicity. Summing over helicities, we obtain 



V. CROSS SECTIONS FOR POLARIZED tf> 



a(qq^^(X)) = 5(s - 4m 2 c ) 




(5.1) 



a(qq — > — 5(s — 4m^) 




(OtCSx)). 



(5.2) 



27ml 
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The cross section for producing a ip with helicity A via the parton process gg — > cc is 
obtained by inserting the matrix elements ( |4.3| ), (|4.5|), and ( |4.13| ) into ( |3.21| ): 



lU8m^ 8 



+ 3(3 -f A °W Po)) + 45(3 "f AO) ^( 3 ^o)) 
mi smi 



27(1 -£ A o) (n i>(ip\\ 

Summing over helicities, this reduces to 

<rto-*) = Ks-4'n>)^l(OiCS )) + "(O*('5 )> + ^<Of( 3 P )} 

Note that the contribution from the matrix element (x ( — ^D 3 )T a i/jV^^o) ip '(— ^D 3 )T a x) 
drops out of the cross section for unpolarized production. If we retain only those matrix 
elements that are of leading order in v 2 , the cross section (|5.3|) reduces to 

a(gg - 4>(\)) - S(s - 4m 2 ) |-^f { (O? (^o)) + 3(3 ~ f A ° W ^o)) } . (5.5) 

288m£ [ J 

The (OjC 1 ^)) term gives V>'s that are unpolarized, while the (Of ( 3 Po)) term gives ^'s with 
helicity —1,0, and +1 in the proportions 3:1:3. 

The production of J/ip from the color-octet parton processes qq — > cc and gg — > cc has 
been studied previously by several groups 0-0. In particular, the spin alignment of the 
ip from these processes has been calculated by Tang and Vantinnen || and by Cho and 
Leibovich |7j. We proceed to compare our results with this previous work. 

Tang and Vantinnen || have used the covariant projection method to calculate the cross 
sections for polarized ip and ip' in fixed-target hadron-hadron collisions from the order-a 2 . 
parton processes. Our result (|5.1| ) for the cross section from qq — » cc is in agreement with 
theirs. Our result ( |5.5| ) for the cross section from gg — > cc agrees with theirs only after 
summing over helicities. We disagree on the helicity dependence of the (0$( 1 So)) term. 
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Their result for this term is proportional to 1 — 5\o, which implies that the ip is produced 
with transverse polarization. This is clearly incorrect, because a pointlike cc pair in a 1 So 
state is rotationally invariant in its rest frame. The t/>'s that are produced by the binding of 
the cc pair must therefore be unpolarized, as in our result ( |5.5| ). 

Cho and Leibovich 0] have used the covariant projection method to calculate the cross 
sections for ip and ip' in pp collisions, including not only the order-a;? parton processes 
but also the order-a^ processes of the form ij — > cc + k. Our results (|5.1|) and ( |5.5|) for 
the cross sections from order-af parton processes agree with theirs. Cho and Leibovich 
determined the short- distance coefficients by matching cross sections for producing color- 
octet cc pairs with vanishing relative momentum and in specific angular momentum states. 
This method gives the correct cross sections for unpolarized production, but it fails in general 
for polarized production. In particular, as we show below, it fails to give the correct short- 
distance coefficients for spin-triplet P-wave matrix elements. In the matrix element (|4.13b|) , 
the indices i and m can be decomposed into contributions from total angular momentum 
J = 0, 1, 2 by using the identity 

feU = UjjJlp fela; l(3\Jh)(Jh\la'; 1/3')) U^U^, (5.6) 

V Jh ) 

where a sum over a, (3, a', and f3' is implied. A similar identity can be used to decompose 
the indices j and n into contributions from total angular momentum J' = 0, 1, 2. The matrix 
element ( [4.13b| ) is thereby expressed as a linear combination of matrix elements of the form 
{J^ (J, h)V^(\)J'(J' , h')) , where the operator J(J,h) creates a pointlike cc pair with total 
angular momentum quantum numbers J and h. The projection operator 7-Wa) projects 
onto states that contain a ip with helicity A plus soft hadrons. If the soft hadrons have total 
angular momentum Jg, then the angular- momentum structure of the matrix element is 

^2(Jh\lX; J S A 5 )(1A; J s X s \J'h') (5.7) 

As 

Rotational symmetry requires that h = h' and that both J and J' lie in the irreducible 
decomposition of 1® Js- It does not require that J = J', so there can be interference terms 
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involving matrix elements with J ^ J'. If we sum (|5.7|) over the helicities A, the orthogonality 
relations for Clebsch-Gordan coefficients imply that J = J'. Thus the interference terms 
cancel upon summing over helicities. The method used by Cho and Leibovich therefore 
gives the correct cross sections for unpolarized production. In the specific case of gg — > cc, 
Cho and Leibovich found that the only amplitudes for the production of spin-triplet P-wave 
cc pairs that are nonzero are gg — > cc(8, 3 Po) and gg — > cc(8, 3 P2, \h\ = 2). Since there is 
at most one nonzero amplitude for any given value of h, there cannot be any interference 
terms. Thus their method gives the correct cross sections for polarized production from 
the order-a^ parton processes. However, for the order-a 3 parton processes ij — > cc + k, 
there are nonvanishing amplitudes with different values of J and the same helicity h, so 
the interference terms are nonzero. These interference terms must be taken into account in 
calculating the term proportional to (Of ( 3 -Po)) in the inclusive cross section for polarized ip. 
The only calculation thus far in which these interference terms have been correctly included 
is a calculation of the gluon fragmentation function for the production of a longitudinally 
polarized ip by Beneke and Rothstein [T^ . 

Fleming and Maksymyk || recently presented a method for calculating cross sections 
for the production of unpolarized heavy quarkonium. To determine the short-distance co- 
efficients of the NRQCD matrix elements, they matched cross sections for producing cc 
pairs with specific color and angular momentum quantum numbers. The cross sections were 
expressed as integrals over the relative 3-momentum q of the cc pair, and they used the non- 
relativistic expansions of the Dirac spinors given in appendix A to expand the integrands as 
Taylor series in q. In general, their method would fail if applied to cross sections for polar- 
ized production, because it does not give the interference terms between parton processes 
that produce cc pairs with different total angular momentum J. Fleming and Maksymyk 
applied their method to the production of unpolarized ip from the parton processes qq —>■ cc 
and gg —>■ cc. Our cross sections agree with theirs after summing over helicities. 
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VI. SUMMARY 



The NRQCD factorization formalism [|l|] is a powerful tool for analyzing the production 
of heavy quarkonium. Cross sections are factored into short-distance coefficients which 
can be calculated using perturbative QCD and long-distance matrix elements which scale 
in a definite way with v. In this paper, we presented a general matching prescription 
for calculating the short- distance coefficients in the inclusive cross sections for quarkonium 
states. Using this matching prescription, the cross sections are obtained in a form that 
applies equally well to any quarkonium state. The specific state enters only in the reduction 
of the NRQCD matrix elements using rotational symmetry and heavy-quark spin symmetry. 
For simplicity, we discussed the reduction of the NRQCD matrix elements only for the case 
of the ip. The generalization to other states will be presented elsewhere ||. 

Our approach has interesting implications for cross sections for producing quarkonium 
states with definite polarization. We showed that cross sections for polarized production 
can involve new matrix elements that do not contribute to cross sections for unpolarized 
production. We also showed that there are interference terms involving the production of 
cc states with different total angular momentum J. These interference terms must be taken 
into account in the cross sections for polarized production from order-af parton processes. 

Our method has a significant advantage over the covariant projection method in that it 
can be straightforwardly generalized to N space dimensions. This is important, because it 
allows the use of dimensional regularization in calculations of radiative corrections. There 
are potential inconsistencies in combining the covariant projection method with dimensional 
regularization for calculations involving orbital angular momentum L = 1 or higher, because 
the projections onto states with definite total angular momentum J are specific to 3 dimen- 
sions. Our method for calculating the short- distance coefficients of NRQCD matrix elements 
involves matching Taylor expansions in the relative momentum, which can be readily calcu- 
lated in N dimensions. After using renormalization to remove the poles in 1/(N — 3) from 
the short-distance coefficients, we can specialize to N = 3 dimensions, and use rotational 
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symmetry and heavy-quark spin symmetry to simplify the matrix elements. Our approach 
thus allows the convenience of dimensional regularization to combined with the full power 
of the NRQCD factorization approach. 

While this paper was being completed, Beneke and Rothstein || presented a paper that 
also points out that there are interference terms involving cc pairs with different total angular 
momentum. They presented a thorough phenomenological analysis of the production in 
fixed-target hadron-hadron collisions of if), if)', and XcJ via order-ct^ parton processes. 
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APPENDIX A: NONRELATIVISTIC EXPANSION OF SPINORS 



In this Appendix, we give the nonrelativistic expansions for the spinors of a heavy quark 
c and antiquark c with arbitrary total 4-momentum P. We assume that the relative 3- 
momentum q of the c in the center-of-momentum (CM) frame of the cc pair is small compared 
to the quark mass m c . The 4-momenta p and p of the c and c can be written 

p = \P + Lq , (Ala) 
p = \P - Lq , (Alb) 

where P is the total 4-momentum and L is a Lorentz boost matrix. ^From the mass-shell 



conditions, p 2 = p 2 = m 2 , we have P ■ Lq = and P 2 = 4E 2 , where E q = \Jm 2 + q 2 . The 
components of the 4-momenta P and Lq in the CM frame of the pair are 



(Lq)" 



= (2E g , 0) , (A2a) 

CM 

= (0, q) . (A2b) 

CM 



When boosted to an arbitrary frame in which the pair has total 3-momentum P, these 
4-momenta are 

P» = (JAB} + P 2 , P) , (A3a) 
(Lq)" = L$q>. (A3b) 

The boost matrix L"-, which has one Lorentz index and one Cartesian index, has components 



1 

pi pj pO pi pi 
2E n P 2 



L\ = —PJ , (A4a) 



The boost tensor L" has many useful properties. Its contraction with the Lorentz vector 
P vanishes: 

P^ 3 = 0. (A5) 
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The contractions of two boost matrices in their Lorentz indices or in their Cartesian indices 
have simple forms: 



-5 ij , 



-<r + 



pppv 



(A6a) 
(A6b) 



There are also simple identities involving contractions of boost matrices with the Levi-Civita 
tensors tp V x P and e ljk : 



, rf tv t\ 

tpuXp L> k 


-. Jjk p 
c J pi 


(A7a) 


rC tv AA 
C/ii/Ap *-i i-L j-T 


— t ijpfc , 


(A7b) 


, rf tv 
tpvXp J- 1 i 1 ^ j 




(A7c) 


, rf pv 
tpvXp -L i-r 


= e lJ,fe L X jL pk , 


(A7d) 


, rC 
tpvXp ^ i 


— e ijk (J=> u L x .L pk -\- P x L p jL uk + P p L U jLxkj , 


(A7e) 


tpvXpP^ 


= e^ k LviL\jL p k, 


(A7f) 



where P' 1 = P^/y/P^. Our sign convention is e *^ = —e 0i jk = e^ k . 

The representation for gamma matrices that is most convenient for carrying out the 
nonrelativistic expansion of a spinor is the Dirac representation: 



7 



-1 



' ^ 



J 



-a 1 



(A8) 



In the CM frame of the cc pair, the spinors for the c and the c are 

1 



u(p) 



v(p) 



CM 



CM 



1 



) 



\ 



-qa-r] 

(E q + m c ) rj ^ 



(A9a) 



(A9b) 



Color and spin quantum numbers on the Dirac spinors and on the 2-component Pauli spinors 
£ and rj are suppressed. When boosted to a frame in which the pair has total 3-momentum 
P, the spinors for the c and c are 
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u(p) = 



v(p) = 



1 



AE q (E P + 2E q )(E q + m c ) 
1 

AE q (E P + 2E q )(E q + m c ) 



(2E q + #> 7o ) 



[2E q +lj> 7o ) 



v q ^ y 

-q- cr r] 

(E q + m c )r] J 



(AlOa) 



(AlOb) 



These spinors are normalized so that uu = —vv = 2m c if the Pauli spinors are normalized 
so that £t£ = ryr] = 1. 

There are 16 independent quantities that can be formed by sandwiching Dirac matrices 
between u{p) and v{p). They are 



u(p)v(p) = -2^(q-a)r), 
u(p)Vv(p) = (W, - P"L",) ( ^ ^Vr? " 



vm r 



E n 



E q (E q + m c ) 



q 3 £}{(\-a)r] 



2 L^ fc e»" H ^ 



«(p)7 m 75^(p) 



P»? V - 2% L M - £*(q x <r) j r], 



(Alia) 
(Allb) 

(Allc) 
(Alld) 
(Alle) 



/From these expressions, it is trivial to carry out the nonrelativistic expansions in powers of 
q. For example, through linear order in q, we have 



u(p)v(p) = -2^(q-a-)r], 

u(p)'f tJ 'v(p) ~ 2m c Ifj £}a^r] , 



( W, - P V L^) - 2 L\ e>"« q l ^V, 



u(ph^v(p) « P M £ f ?7 - 2i L^- ^(q x cr) j r}, 



(A12a) 
(A12b) 
(A12c) 
(A12d) 
(A12e) 
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APPENDIX B: PRODUCTION MATRIX ELEMENTS 



In this Appendix, we define some of the standard NRQCD production matrix elements 
that were introduced in Ref. M. In order to establish the notation for the fields, we give the 
lagrangian for minimal NRQCD: 

c = if (w, + ^) * + x< (iD, - i£) X + C QCD , (Bl) 

where tp is the Pauli spinor field that annihilates a heavy quark, x is the Pauli spinor field 
that creates a heavy antiquark, and Cqcd is the usual QCD lagrangian for the gluons and 
the light quarks and antiquarks. 

The production operators in Ref. jl| were defined using a projection operator a^ H an = 
p m?a) defined by 

-Pnfx) = y £\H(P = 0,\) + S)(H(P = 0,\) + S\, ( B2 ) 

s 

where the states in the sum have the standard nonrelativistic normalization. For example, 
the normalization of the quarkonium state is 

(H(P', \')\H(P, A)) = (2vr) 3 5 3 (P - P') 5 XX , (B3) 

Thus the projection operator (|B!|) has energy dimension -3. 
The production operators of dimension 6 are 

O? CS ) = xV Kfx) <A f X, (B4a) 

Of ( 3 ^) = X VV V$fa V>V X , (B4b) 

OffSo) = x ] T a ^ 7^f A) ^T a x , (B4c) 

Oi ( 3 ^) = X VT> pNRj ^t a i T a x _ (B4d) 

Some of the production operators of dimension 8 are 

OfC^) = X f (-iB> V%f x) ^(-0%, (B5a) 



1 



Of ( 3 Po) = 77 X f H D • o-)^ P£f A) ^ f (-f D • (B5b) 
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Of ( l P 1 ) = x t ( _0* )T> P NR } ^t ( _0 !)rXj (B5c ) 
O^'Po) = I X f (-|D ■ (r)T a i/j V%f x) ^ f (-fD ■ <r)T a X , (B5d) 

where = ^Dx — (DipYx- 

Note that the projection operator defined by (|B2|) differs from the corresponding projec- 
tion operator ( |2.4| ), which is defined using states with the standard relativistic normalization. 
At leading order in v 2 , these operators differ simply by an overall factor: Vh(\) ~ ^ m cP^\y 
Beyond leading order, the relation between the two projectors is complicated because the 
normalization differs for each term in the sum over soft states S. For physical quantities, like 
the cross section in ( |3.1| ), the difference between the projection operators is compensated by 
the short- distance coefficients: 

E Cmn {X^'U VH(X) V f £nX> = E (X^'U Kfx) ■ (B6) 

mn ran 

We now list the relations between scalar matrix elements defined with the projection 
operator ( |2.4j) and the matrix elements of the operators defined above. At leading order in 



a s and to leading order in v 2 , we have 

(xH Vh ( x) tfx) = 4 m c O?CS ), (B7a) 

(xVV Ph W V>V X > = 4 m c Of ( 3 ^), (B7b) 

( x t(_i£i)^ (A) ^t(_i£i )x ) = 4m c Of( 1 P 1 ), (B7c) 

(^(-fD • a)i> Pfjfa V f (-|D • <r)x> = 12 m c Of ( 3 P ). (B7d) 

The corresponding color-octet matrix elements are 

( X f T> 7> H(A) ^T a X ) = 4 m c Of (^a), (B8a) 

<*VT> p ff(A) V>VT a X > = 4 m c Of ( 3 Sx), (B8b) 

( x t(_i^)T> Ph ( a) ^(-fS J )T a X ) = 4m c Of( 1 P 1 ), (B8c) 

(X f (-1D • <r)T> P*f A) ^(-ID • o-)T a X ) = 12 m c Of ( 3 P ). (B8d) 
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